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Abstract 

We report evidence for direct CP violation in the decay B° K'^Tr~ with 253 fb~^ of data col- 
lected with the Belle detector at the KEKB e~^e~ collider. Using 275 million BB pairs we observe a 
B —>■ K^TT^ signal with 2 140 ±53 events. The measured CP violating asymmetry is Acp{K~^'^~) = 
—0.101 zb 0.025(stat) ± 0.005(syst), corresponding to a significance of 3.9a" including systematics. 
We also search for CP violation in the decays S+ iC+7r° and 5+ tt+tt^. The measured 
CP violating asymmetries are AcpiK^-^^) = 0.04 ± 0.05(stat) ± 0.02(syst) and ^cp(7r+7r°) = 
-0.02 ± O.lO(stat) ± O.Ol(syst), corresponding to the intervals -0.05 < AcpiK'^T^^) < 0.13 and 
-0.18 < ^cp(vr+7r°) < 0.14 at 90% confidence level. 

PACS numbers: ll.30.Er, 12.15.Hh, 13.25.Hw, 14.40.Nd 
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In the Standard Model (SM), CP violation arises via the interference of at least two 
diagrams with comparable amplitudes but different CP conserving and violating phases. 
Mixing induced CP violation in the B sector has been established in b ^ ccs transitions 

0]- In the SM, direct CP violation is also expected to be sizable in the B meson system 
3|. The first experimental evidence for direct CP violation in B mesons was shown by Belle 
br the decay mode B^ tt+tt" This result suggests large interference between tree 
and penguin diagrams and the existence of final state interactions j^. Recently, both Belle 
PI and BaBar [3] have reported searches for direct CP violation in another decay mode 
B^ K^Ti^ , where direct CP violation is also expected. 

The CP violating partial rate asymmetry is measured as: 

^ N(E ^ 7) - NjB ^ /) 
N(B-.f) + N{B^fy 

where N{B —>■ f) is the yield for the B — > Kh/'kit decay and N{B — > /) denotes that of 
the charge-conjugate mode. Theoretical predictions with different approaches suggest that 
Acp{K^T^^) could be either positive or negative [1]. Although there are large uncertain- 
ties related to hadronic effects in the theoretical predictions, results for Acp{K^t^~) and 
Acp{K^TT^) are expected to have the same sign and be comparable in magnitude p]. In this 
Letter, we report Acp measurements for B^ — > /T^vr", B^ K^ti^ and B^ — >■ 7r"'"7r° using 
275 million BB pairs collected with the Belle detector at the KEKB e^e~ asymmetric-energy 
(3.5 on 8 GeV) collider operating at the T(4S') resonance. 

The Belle detector is a large-solid-angle magnetic spectrometer that consists of a silicon 
vertex detector (SVD), a 50-layer central drift chamber (CDC), an array of aerogel threshold 
Cherenkov counters (ACC), a barrel-like arrangement of time-of-flight scintillation counters 
(TOF), and an electromagnetic calorimeter (ECL) comprised of CsI(Tl) crystals located 
inside a superconducting solenoid coil that provides a 1.5 T magnetic field. An iron flux- 
return located outside of the coil is instrumented to detect mesons and to identify muons 
(KLM). The detector is described in detail elsewhere Two different inner detector 

configurations were used. For the first sample of 152 million BB pairs (Set I), a 2.0 cm 
radius beampipe and a 3-layer silicon vertex detector were used; for the latter 123 million 
BB pairs (Set II), a 1.5 cm radius beampipe, a 4-layer silicon detector and a small-cell inner 
drift chamber were used ^ 



The B candidate selection is the same as described in Ref. Charged tracks are 

required to originate from the interaction point (IP). Charged kaons and pions are identified 
using dE/dx information and Cherenkov light yields in the ACC. The dE/dx and ACC 
information are combined to form a K-ir likelihood ratio, TZ{Ktt) = £k/{^k + ^n), where 
jC-k (-Ctt) is the likelihood of kaon (pion). Charged tracks with IZi^K-w) > 0.6 are regarded 
as kaons and tracks with TZ{Ktx) < 0.4 as pions. Furthermore, charged tracks that are 
positively identified as electrons are rejected. The electron identification uses the information 
composed of E /p and dE/dx, shower shape, track matching x^, and ACC light yields. The 
K/ir identification efficiencies and misidentification rates are determined from a sample of 
kinematically identified D*~^ D^n^^D^ K~tt~^ decays, where the kaon and pion from 
the D decay are selected in the same kinematic region as in the B^ K^tt" decay. Table E] 
shows the results. The detection efficiency for K~7t'^ is found to be 1.0% greater than that 
for K~^7r~] this small difference is corrected for in the Acp measurement. 

Candidate vr^ mesons are reconstructed by combining two photons with invariant mass be- 
tween 115 MeV/c^ and 152 MeV/c^, which corresponds to ±2.5 standard deviations around 
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TABLE I: Performance oi K — it identification measured using D*+ D^ir^, K tt^ decays. 



Set I Set II 

Eff. (%) Fake rate (%) Eff. (%) Fake rate yvo 

if+ 83.74 ±0.18 5.10 ±0.12 82.41 ± 0.20 6.57 ±0.15 

ET" 84.73 ±0.18 5.69 ±0.12 83.26 ± 0.20 7.14 ±0.15 

7r+ 91.25 ±0.15 10.74 ±0.15 89.48 ±0.18 11.82 ±0.17 

TT- 90.54 ±0.16 10.09 ±0.15 88.56 ±0.19 11.57 ±0.17 



the nominal 7r° mass. Each photon is required to have a minimum energy of 50 MeV in 
the barrel region (32° < d., < 129°) or 100 MeV in the end-cap region (17° < 6^ < 32° 
or 129° < 6^ < 150°), where 6.^ denotes the polar angle of the photon with respect to the 
beam line. To further reduce the combinatorial background, 7r° candidates with small decay 
angles (cos^* > 0.95) are rejected, where 9* is the angle between the vr" boost direction in 
the laboratory frame and its 7 daughters in the 7r° rest frame. 

Two variables are used to identify B candidates: the beam-constrained mass, Mbc = 

\/-^beam ~ ' ^"^^ energy difference, = — -^beam' where -Ebeam is the beam 
energy and E]^ and p*^ are the reconstructed energy and momentum of the B candidates in 
the center-of-mass (CM) frame. Events with Mbc > 5.20 GeV/c^ and -0.3 GeV < < 
0.5 GeV are selected for the final analysis. 

The dominant background is from e~^e~ qq {q = u, d, s, c) continuum events. To 
distinguish the signal from the jet-like continuum background, event topology variables and 
B flavor tagg ing information are employed. We combine a set of modified Fox- Wolfram 
moments [13| into a Fisher discriminant. The probability density function (PDF) for this 
discriminant, and that for the cosine of the angle between the B flight direction and the 
z axis, are obtained using signal and continuum Monte Carlo (MC) events. These two 
variables are then combined to form a likelihood ratio 71 = Cs/{Cs + Cqq), where Cs{qq) is 
the product of signal {qq) probability densities. Additional background discrimination is 
provided by B flavor tagging. The standard Belle flavor tagging algorithm gives two 
outputs: a discrete variable indicating the flavor of the tagging B and a MC-determined 
dilution factor r, which ranges from zero for no flavor information to unity for unambiguous 
flavor assignment. An event with a high value of r ( typically containing a high- momentum 
lepton) is more likely to be a BB event so a looser TZ requirement can be applied. We divide 
the data into r > 0.5 and r < 0.5 regions. The continuum background is reduced by applying 
a selection requirement on TZ for events in each r region of Set I and Set II according to 
the figure of merit defined as N^^^/ Ns^^ + Ngq^, where N^^p denotes the expected signal 
yields based on MC simulation and our previous branching fraction measurements and 
Nq^^ denotes the expected qq yields from sideband data (Mbc < 5.26 GeV/c^). A typical 
requirement suppresses 92-99% of the continuum background while retaining 48-67% of the 
signal. 

Backgrounds from T(45') —>■ BB events are investigated using a large MC sample. Af- 
ter the TZ requirement, we find a small charmless three-body background at low AE, and 
reflections from B^ — > 7c'^7i~ decays due to K-tt misidentification. 

The signal yields are extracted by applying unbinned two dimensional maximum likeli- 
hood (ML) fits to the (Mbc and AE) distributions of the B and B samples. The likelihood 
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for each mode is defined as 



L = exp (- J2 Ns,k,j) X n(E Ns,k,jVs,k,j,i) (2) 

s,k,j i s,k,j 

'Ps,k,j,i = ^[l - qiAcpj]Ps,k,j{Mkci,AEi), (3) 

wliere s indicates Set I or Set II, k distinguishes events in the r < 0.5 or r > 0.5 regions, i 
is the identifier of the i-th event, P(Mbc, AE) is the two-dimensional PDF of Mbc and AE, 
q indicates the B meson flavor, B{q = +1) or B{q = —1), Nj is the number of events for 
the category j, which corresponds to either signal, qq continuum, a reflection due to K-tt 
misidentification, or background from other charmless three-body B decays. 

The yields and asymmetries for the signal and backgrounds are allowed to float in all 
modes. Since the K'^ir^ and ir'^n^ reflections are difficult to distinguish with AE and Mbc 
we fit these two modes simultaneously with a fixed reflection-to-signal ratio based on the 
measured K-n identification efficiencies and fake rates. All the signal PDFs (P(Mbc, Ai?)) 
are obtained using MC simulations based on the Set I and Set II detector configurations. The 
same signal PDFs are used for events in the two different r regions. No strong correlations 
between Mbc and AE are found for the B K^tt~ signal. Therefore, its PDF is modeled 
by a product of a single Gaussian for Mbc and a double Gaussian for AE. For the modes 
with neutral pions in the final state, there are correlations between Mbc and AE in the tails 
of the signals; hence, their PDFs are described by smoothed two-dimensional histograms. 
Discrepancies between the signal peak positions in data and MC are calibrated using B^ — »• 
D^n'^ decays, where the — > K~^n~7r^ sub-decay is used for the modes with a tt^ meson 
while K^Ti^ is used for the K^7t~ mode. The MC-predicted AE resolutions are 

verified using the invariant mass distributions of high momentum D mesons. The decay 
mode — > K~^7i~ is used for B^ — > K~^7i~, and — >■ K^ir'n^ for the modes with a 7r° in 
the final state. The parameters that describe the shapes of the PDFs are fixed in all of the 
fits. 

The continuum background in AE is described by a first or second order polynomial while 
the Mbc distribution is parameterized by an ARGUS function f{x) = xx/l — exp [—^(1 — 
x^)], where x is Mbc divided by half of the total center of mass energy [l^. The continuum 
PDF is the product of an Argus function and a polynomial, where ^ and the coefficients 
of the polynomial are free parameters. These free parameters are r-dependent. A large 
MC sample is used to investigate the background from charmless B decays and a smoothed 
two-dimensional histogram is taken as the PDF. The functional forms of the PDFs are the 
same for the B and B samples. 

The efficiency of particle identification is slightly different for positively and negatively 
charged particles; consequently the raw number of asymmetry in Eq. 3 no longer gives Acp 
correctly and must be corrected. For the K'^tt^ and ir'^n^ modes, this raw asymmetry can 
be expressed as: 

- 1 + A.Acp' 

where Acp is the true partial rate asymmetry and the efficiency asymmetry Ae is the 
efficiency difference between K~{7r~^) and K~^{Tr~) divided by the sum of their efficiencies. 
The situation is more complicated for the K~^tt~ mode because, in addition to the bias 
due to the efficiency difference between K~n^ and K~^7t~, a K~7i~^ signal event can be 
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AE (GeV) AE (GeV) 



FIG. 1: Mbc (top) and AE (bottom) distributions for B K-Tr+ (left) and 5° K+tt' (right) 
candidates. The histograms represent the data, while the curves represent the various components 
from the fit: signal (dashed), continuum (dotted), three-body B decays (dash-dotted), background 
from mis-identification (hatched), and sum of all components (solid). 

misidentified as a K~^7r~ candidate and dilute Acp- The efficiency asymmetry results in an 
Acp bias of +0.01 while the small dilution factor due to double misidentification reduces 
•Acp by a factor of 0.99. These effects are included in the raw asymmetry correction and 
their errors are included in the systematic uncertainty. 

Table HTl gives the signal yields and Acp values for each mode. The asymmetries for the 
background components are consistent with zero within errors. Projections of the fits are 
shown in Figs. ^El The systematic errors from fitting are estimated from the deviations 
in Acp after varying each parameter of the signal PDFs by one standard deviation. The 
uncertainty in modeling the three-body background is studied by excluding the low AE 
region (< —0.12 GeV) and repeating the fit. Systematic uncertainties due to particle iden- 
tification are estimated by checking the fit after varying the K/tt efficiencies and fake rates 
by one standard deviation. At each step, the deviation in Acp is added in quadrature to 
provide the systematic errors, which are less than 0.01 for all modes. A possible bias from 
the fitting procedure is checked in MC and a bias due to the 71 cut is investigated using the 
D^7i~^ samples. No significant bias is observed. The systematic uncertainties due to 
the detector bias are obtained using the fit results for the continuum background listed in 
Table ITTl The final systematic errors are then obtained by quadratically summing the errors 
due to the detector bias and the fitting systematics. 

The partial rate asymmetry Acp{K~^tt~) is found to be —0.101 ± 0.025 ± 0.005, which 
is 3.9a from zero. The significance calculation includes the effects of systematic uncer- 
tainties. Our result is consistent with the value reported by BaBar, Acp{K~^it~) = 
—0.133 ± 0.030 ± 0.009 0]. The combined experimental result has a significance greater 
than 5(7, indicating that direct CP violation in the B meson system is established. Our 
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FIG. 2: Mbc (top) and AE (bottom) distributions for B' R-tt^ (left) and B+ K+tt^ (right) 
candidates. The curves are described in the caption of Fig. ^ 

TABLE II: Fitted signal yields, Acp results and background asymmetries for individual modes. 



Mode 


Signal Yield 


Acp 


Bkg Acp 




2140 ± 53 


-0.101 ±0.025 ±0.005 


-0.001 ± 0.005 




728 ± 34 


0.04 ±0.05 ±0.02 


-0.02 ±0.01 




315 ±29 


-0.02 ±0.10 ±0.01 


-0.01 ±0.01 



measurement of Acp{K^tt^) is consistent with no asymmetry; the central value is 2Aa 
away from AcpiK^ir'). If this result is confirmed with higher statistics, the difference may 
be due to the contribution of the electroweak penguin diagram or other mechanisms 
No evidence of direct CP violation is observed in the decay -k^-k^ . We set 90% CL. 

intervals -0.05 < Acp{,K+'k^) < 0.13 and -0.18 < Acp{,-n^-n^) < 0.14. 
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FIG. 3: Mbc (top) and AE (bottom) distributions for B' tt'tt^ (left) and 
candidates. The curves are described in the caption of Fig. ^ 
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